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摘要 ”相对论 性 重 离子 对 撞 的 目的 之 一 是 寻找 量子 色 动 力学 (Quantum Chromodynamics , QCD) ) 临 界 终点 
(Critical End Point, CEP) ,实验 中 测 到 的 净 质 子 数 扰动 呈现 出 非 单 调 的 行为 ,这 暗示 了 CEP 的 存在 。 本 文 使 用 
了 3 味 PNJL(Polyakov-loop Nambu-Jona-Lasinio) 模 型 , 沿 着 从 实验 数据 中 拟 合 的 化 学 冻结 线 ,计算 了 重子 数 扰 
动 累积 量 之 比 CyC, 随 对 撞 能 量 的 变化 。 结 果 发 现 , 基 于 平衡 态 假设 ,在 对 撞 能 量 7.7~200 GeV 范围 内 , 随 着 对 
撞 能 量 的 降低 ,CYC, 先 缓慢 下 降 后 上 升 , 这 与 实验 数据 一 致 。 这 也 暗示 平衡 态 假设 可 以 用 于 探索 重 离子 对 撞 后 
系统 的 演化 行为 ,揭示 了 相 变 线 跟 化 学 冻结 线 的 关系 对 观测 量 的 影响 非常 重要 。 

关键 词 QCD 相 图 ， 临 界 终 点 ， 重 离子 对 撞 ， 重 子 数 扰动 
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QCD critical end point and baryon number fluctuation 


XU Kun HUANG Mei 


(School of Nuclear Science and Technology, University of Chinese Academy of Sciences, Beijing 101400, China) 


Abstract One of the main goals of relativistic heavy-ion collision (HIC) is to search for the critical end point (CEP) 
of quantum chromodynamics (QCD), and distribution of the net-proton number from experimental measurements 
shows non-monotonic behavior, which indicates the existence of a CEP. The purpose of this work is to investigate the 
relationship between the net-proton number fluctuation and collision energy, and to explain the experimentally 
measured behavior. This study investigates the three-flavor Polyakov-loop Nambu-Jona-Lasinio (PNJL) model, 
which contains quark degrees from the NJL (Nambu-Jona-Lasinio) model and effective gluon contributions from 
Polyakov-loop, based on the equilibrium assumption and mean-field approximation. In addition, we study the phase 
diagram and C,/C, of baryon number fluctuation as a function of collision energy along the freeze-out lines fitted 
from experimental data. With an appropriate form of freeze-out line, the collision energy decreases in the region of 
7.7~200 GeV, and C/C, decreases slightly then increases, which is in agreement with the experimental data. 
Additionally, these results indicate that the equilibrium assumption is appropriate for the exploration of the system 
evolution after HIC, and the relationship between the freeze-out and phase transition lines is highly sensitive for 


observables. 
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构 是 近 些 年 的 研究 热点 ,在 极端 条 件 下 ,QCD 物质 
会 呈现 出 完全 不 同 于 真空 下 的 形态 。 随 着 温度 的 升 
高 ,QCD 物质 会 从 强 子 气态 变 为 夸克 胶 子 等 离子 体 
(Quark-Gluon Plasma, QGP) ,自由 度 也 从 强 子 转 变 
为 夸克 和 胶 子 , 格 点 QCD (Lattice QCD ) 发 现在 小 化 
学 势 时 手 征 相 变 与 退 禁 闭 相 变 均 为 光滑 过 渡 
(Crossover)"。 另 一 方面 ,从 对 称 性 分 析 以 及 手 征 
有 效 模 型 计算 ,人 们 普遍 认为 在 高 化 学 势 处 发 生 一 
级 手 征 相 变 。 因 此 ,根据 朗 道 相 变 理论 ,在 有 限 温 度 
有 限 化 学 势 处 则 存在 一 级 相 变 终点 , 即 QCD 临界 终 
点 (Critical End Point, CEP)". 

高 温 高 密 的 极端 环境 与 早期 宇宙 演化 及 天 体 物 
理 有 很 强 的 相关 性 ,高 密 环境 存在 于 致密 星 内 部 , 比 
如 中 子 星 。 强 相互 作用 物质 在 有 限 温度 有 限 密度 时 
的 状态 方程 ,直接 影响 了 中 子 星 质量 -半径 关系 。 在 
实验 室 中 通过 重 离 子 对 撞 (Heavy Ion Collision, 
HIC) 产 生 高 温 高 密 的 环境 ,是 研究 强 相 互 作用 物质 
在 极端 条 件 下 性 质 的 重要 场所 ,比如 研究 轻 核 性 
BOS) it AS DE UT as E SB HE 
Jg 077 WE 55 Ex e PO, DA ASIC ES APER TJ 
QCD fH Fl UL QCD 临界 点 。 

在 重 离子 对 撞 中 ,假设 化 学 冻结 线 与 相 边 界 足 
够 近 , 那 么 相 图 信息 则 可 以 通过 守恒 荷 的 扰动 呈现 
出 来 。 据 预测 ,这 些 守恒 荷 对 关联 长 度 非常 敏感 。 
因此 在 重 离子 对 撞 中 ,通过 分 析 逐 事件 的 重子 数 、 电 
荷 以 及 奇异 数 的 扰动 ,尤其 是 它们 的 累积 量 ,能 够 揭 
示 RHIC (Relativistic Heavy-ion Collision) 以 及 LHC 
(Large Hadron Collider) 的 重 离子 对 撞 中 产生 的 物 
SA ATI AE. SET EE FIE tat Et IN 79— X 
对 撞 事 件 的 测量 结果 ,在 大 量 记录 事件 中 ,此 守恒 蓓 
的 累积 量 (cumulant) 定 义 如 下 : 


N 


C= (N) 
C, - (GN) 
C, = ((ÓN y) ( ) 


C,=((6N)*) - 3((N y) 
其 中 :5N = N -《N》, 这 些 累积 量 与 不 同 矩 直接 的 关 
C, C, 


M Ce C,,8 (C yr" — (C y (2) 
2 2 


其 中 :Mo 分 别 为 对 应 守恒 和 荷 的 平均 值 
(Mean) , 77 4 (Variance) ~ {iri S£ (Skewness) VA Ae we BE 
(Kurtosis)。 在 实验 中 ,累积 量 是 与 系统 体积 成 正比 
的 ,为 了 消除 体积 带 来 的 影响 ,人 们 往往 使 用 累积 量 
的 比值 ,例如 : 


C C 
2 = 2 = 3 
o’ IM C. So "x KG C, 


(3) 


有 了 以 上 的 定义 ,就 可 以 在 给 定 的 对 撞 能 量 
Vsw 下 ,通过 测量 逐 事件 的 净 质 子 分 布 计算 各 个 累 
积 量 及 其 比值 。 

RHIC 的 能 量 束 流 扫 描 计划 (Beam Energy 
Scan,BES) 通 过 金 - 金 对 撞 ,测量 了 净 质 子 数 、 净 电 
荷 数 以 及 净 及 介子 数 分 布 的 累积 量 , 对 撞 能 量 Vs 
为 7.7 GeV 11.5 GeV.14.5 GeV 、19.6 GeV.27 GeV. 
39 GeV、54.4 GeV.62.4 GeV 以 及 200 GeV, 其 结果 
见 文献 [36-39]。 实 验 结果 显示 ,在 金 - 金 中 心 对 撞 
中 , 净 质 子 数 分 布 的 C,/C, 出 现 了 随 对 撞 能 量 的 非 
单调 的 行为 :在 对 撞 能 量 为 200 GeV HY ,C,/C, AA 
1, 随 着 对 撞 能 量 的 下 降 , C/C, 先 在 对 撞 能 量 为 
20 GeV 时 降低 到 0.1, 然 后 又 在 对 撞 能 量 为 7.7 GeV 
时 迅速 升 高 到 3.5 左 右 。STAR 最 新 的 3 GeV AEH 
实验 结果 以 及 HADES 的 2.4 GeV 对 撞 都 显示 了 
Cs/C; 为 负 值 中 。 从 以 上 的 实验 数据 中 观测 到 重子 
数 扰动 的 震荡 行为 ,我 们 需要 知道 可 以 从 中 提取 关 
F CEP 的 何 种 信息 以 及 CEP 位 于 相 图 中 的 何 处 。 
一 方面 ,应 该 从 沿 着 化 学 冻结 线 的 测量 结果 中 提取 
KF QCD 相 变 的 有 效 信息 ; 另 一 方面 ,也 应 该 小 心 
地 探索 QCD 相 变 如 何 影响 诸如 守恒 量 的 累积 量 等 
性 质 。 图 1 展示 了 在 3 GeV 对 撞 能 量 以 及 HADES 
在 2.4 GeV 对 撞 能 量 下 的 结果 "”。 


STAR (0%~5%) 


Ratio C/C, 


5 10 50 100 
Collision Energy \ Syy / GeV 


1 STAR 测量 的 金 - 金 对 撞 中 , 净 质 子 数 的 CyC, 

随 对 撞 能 量 的 变化 

Fig.1 Variation of the ratios of cumulants, C,/C,, with 
collision energy for net-proton measured by STAR 


要 理解 实验 中 出 现 的 非 单调 的 行为 ,就 需要 从 
理论 上 研究 中 质子 数 扰动 随 对 撞 能 量 的 行为 。 对 于 
平衡 态 系统 ,一般 在 理论 中 可 以 定义 热力 学 势 , 形 
式 为 : 


Q- En (ZU, Tous) A) 
RE: TAME: VARIAR: ZIRA R: 
pos 分 别 为 重子 数 化 学 势 . 电 荷 数 化 学 势 以 及 
奇异 数 化 学 势 。 有 了 热力 学 势 , 系 统 的 压强 也 就 自 
然 知道 了 :P =- Q. SF TE far TIAUS BE (susceptibility) 
定义 为 无 量 纲 压强 对 相应 化 学 势 的 导数 : 
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"TU 
Kaas 
式 中 :fi T uu/T,q- B,Q,S。 这 些 守 恒 荷 分布 的 累 
积 量 与 对 应 敏感 度 的 关系 为 : 
9 Ig (Z CV, T, un, Ha Hs )) = 

oj ojo 

VI^ ys" (T, Ue Hos Ms) (6) 

通过 这 个 关系 ,就 可 以 在 理论 上 计算 重子 数 扰 
动 随 着 温度 、 化 学 势 以 及 对 撞 能 量 的 变化 ,并 与 实验 
结果 进行 比较 ,进而 提取 关于 QCD 相 变 的 信息 。 

在 不 同 的 模型 中 ,基本 自由 度 及 其 相互 作用 的 
形式 不 同 , 配 分 函数 也 不 同 , 因 此 得 到 的 热力 学 性 质 
以 及 相 变 也 就 不 同 。 人 们 使 用 了 不 同 的 模型 ,在 理 
论 中 对 QCD 临界 点 进行 广泛 地 分 析 过 , 比如 
Nambu-Jona-LasinioCNJL) 模 型 ,使 用 Polyakov-loop 
改进 的 NJL 模型 (PNJL) 、 线 性 c 模 型 、 夸 克 - 介 子 
(Quark-Meson , QM) 模型 、Polyakov-loop 改进 的 夸 
克 - 介 子 模型 .Dyson - Schwinger 7; f£ (DSE) .@& E 
QCD 模型” 。 另 外 ,利用 流体 模型 分 析 重 子 数 扰 
动 以 及 重 离子 对 撞 中 其 他 现象 也 是 目前 主流 的 方 


BQS — 
C 


nmi 2 m; 
Hn = B8 (TT) KT) 
cosh (B, is + Qio + Sis) 
m, m, 
NT = 25 ( T YK, ( T ) X 
sinh (B; ûs + Qio + Sifts) 


因此 ,考虑 质子 ,也 就 是 B = 1, 重 子 数 敏感 度 的 
比率 为 : 


(8) 


(9) 


Cam Keven ] 

C2 B ? 
even X 

C? B 

coe edi (10) 
odd  Xodd 

Cd - 


B 
n = even = coth ( /T )|, — 
Cha Xaa P EL S 


本 文 关 注 的 四 阶 与 二 阶 累积 量 的 比值 则 为 : 


C4 

cM 1 (11) 
也 就 是 说 ,在 低温 低 密 的 强 子 相 , C4/C, = 1, 这 

可 以 作为 强 子 相 的 信号 。 而 随 着 温度 的 升 高 到 相 变 

度 ( 约 150 MeV) , 强 子 逐渐 “融化 ”, 系 统 的 主要 自 

昌 度 变 为 夸克 和 胶 子 ,因此 C,/C, 也 会 逐渐 偏离 1， 


as 
Kg0Og 一 


法 中。 然而 ,不 同 的 模型 ,甚至 是 相同 的 模型 使 用 
不 同 的 参数 ,给 出 的 CEP 的 位 置 都 不 同 中 。 因 此 ， 
结合 理论 与 实验 的 结果 探索 CEP 是 否 存在 甚至 确 
定 CEP 的 位 置 是 必要 的 。 


1 ” 强 子 共振 气 模型 


在 低温 低 密 的 时 候 , 强 子 是 强 相互 作用 物质 的 
主要 自由 度 , 随 着 温度 的 升 高 ,自由 度 就 由 强 子 变 为 
夸克 和 胶 子 。 在 强 子 共振 气 (Hadron Resonance 
Gas,HRG) 模 型 中 ,无 相互 作用 的 强 子 以 及 共振 态 
是 基本 自由 度 。 在 玻 尔 兹 曼 近 似 下 ,系统 的 压 
gm Jet. 


P 1 
T^ 一 yp (Z(V, T, le, Ho Ls )) 
= m, m, (7) 
BEET KT) 


cosh (Bifis + Qila + Sifts) 

式 中 :m, 是 强 子 i 的 质量 ;g, 是 其 简 并 度 ;K, 是 修正 贝 
塞 尔 函 数 。 求 和 遍历 所 有 的 强 子 以 及 共振 态 ,包括 
对 应 的 反 粒 子 。 由 HRG 得 到 的 结果 一 般 可 以 作为 
在 重 离子 对 撞 中 寻找 QCD 相 变 以 及 临界 点 的 基线 。 
对 于 净重 子 数 扰动 来 说 ,累积 量 的 比率 的 形式 非常 
简单 ,在 玻 尔 兹 曼 近似 下 ,重子 数 敏感 度 的 形式 可 以 
写 为 : 


这 也 是 发 生 相 变 的 判 据 之 一 。 
2 格 点 QCD(Lattice QCD) 


格 点 QCD 是 通过 没有 任何 假设 的 第 一 性 原理 
一 一 即 量 子 色 动力 学 进行 计算 的 非 微 扰 方 法 ,完全 
可 以 用 于 计算 处 于 热平衡 系统 的 强 相互 作用 体系 的 
热力 学 性 质 。 由 于 “符号 问题 ”, 格 点 QCD 无 法 直接 
计算 有 限 化 学 势 的 区 域 ,但 是 通过 泰勒 展开 技巧 ,可 
以 近似 计算 小 化 学 势 的 区 域 。 具 体 来 说 ,有 限 化 学 
势 人 时 的 压强 可 以 展开 为 "”: 


POS) PU) = dry By tt 
1 x (T) (y+ Xs (7) (Ey 


12 y?(T) T 360*?(T) T 

O (ug ) (12) 

FHF QCD 的 对 称 性 ,奇数 项 贡献 为 零 ,只 有 偶 
数 项 有 不 为 零 。 可 以 看 到 ,展开 系数 与 重子 数 的 敏 
RE y y Hy S VAR y Ey TES 

同样 的 原理 ,也 可 以 利用 重子 数 敏感 度 计 算 有 
限 化 学 势 下 的 值 : 


g jid Eu 
ge Hoo — qug 
k-ol** 


2 EAN T Ti jx QCD? R Se OF 
a3 lye 随 温 度 的 变化 以 及 压强 在 不 同 的 化 学 势 下 随 
温度 的 变化 。 可 以 看 到 ,在 低温 下 与 HRG 的 结果 符 
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thy 


160 180 200 220 240 260 280 
T/MeV 


2 BNL-Bielefeld-CCNU 通过 格 点 QCD 计算 的 在 零 化 学 
35 F yt? 随 温 度 的 变化 
Fig.2 Lattice QCD result for 2/7? as a function of 
temperature at zero baryon chemical potential from BNL- 
Bielefeld-CCNU collaboration 


Ao, EATEN st iy HALF HRG 模型 的 结 
Rot? = 1。 如 果 要 计算 大 化 学 势 的 区 域 , 则 需 
要 计算 到 更 高 阶 的 结果 ,这 就 需要 更 多 的 计算 ,花费 
更 多 的 时 间 。 目 前 格 点 QCD 的 计算 通过 泰勒 展开 


还 不 能 推广 到 足够 高 的 化 学 势 , 在 能够 应 用 泰勒 展 
开 的 区 域 目前 没有 发 现 CEP 的 存在 。 


3 3 味 Polyakov-loop Nambu-Jona-Lasinio 
模型 


为 了 能 够 研究 有 限 温 度 有 限 化 学 下 的 相 图 ,使 
用 合适 的 有 效 模型 进行 计算 是 非常 必要 的 。NJL 
(Nambu-Jona-Lasinio ) 模 型 经 常 被 用 来 研究 QCD 相 
图 ,但 是 NJL 模 型 给 出 的 相 变 温度 比 格 点 QCD 计算 
结果 高 ,因此 ,选用 3 味 的 Polyakov-loop 改进 的 NJL 
模型 ,这 个 模型 给 出 的 手 征 相 变温 度 跟 格 点 结果 相 
当 , 因 此 被 叫做 实际 的 NJL 模型 (简写 PNJL, Bl 
realistic PNJL) ,TPNJL 模型 考虑 了 8 夸克 相互 作用 。 
这 里 直接 写 出 有 限 温度 、 有 限 化 学 势 下 的 平衡 态 系 
SNA 56340057. 


2 


Q-g. 2 - Er o 0,0, + 3% | > 十 382 D9} - z 2s 
d? P "a (Er- uj)iT | - (E, - nj) 3(E,- nj)! 
ar | ^e +3(@ + de Je . ma 
) p (E,* uj)iT | (Ert ym e 1E + yu) 
ar | rn + x(2 + oe Je "|+ 
ed 
式 中 :第 一 行为 真空 部 分 的 能 量 ; 第 二 行 与 第 三 行 则 UU xin J(®,®) | (16) 
为 夸克 在 有 限 温 度 有 限 化 学 势 下 的 贡献 ;第 四 行为 rT T 
Polyakov-loop 的 贡献 。g,, gm gl, 8 分 别 为 夸克 相互 iem) " " 
作用 的 耦合 常数 :or = (Ty) AS v6 凝聚 ;f= u, d Xt T 5 BP- 6 (d? d») 4 (8P (17) 


应 两 味 轻 夸 克 而 f=s 对 应 奇异 夸克 。Ej= 


Jp + My 为 夸克 的 色散 关系 ,其 中 动力 学 质量 MM 


的 形式 为 : 


M,- m, - 2g,0,+ B oy, 10f+2 7 28,0,» 0} - 4g,0; 
= 


(15) 
如 果 o = co， 那么 co = oo = 0,。 由 于 NIL 
模型 部 分 为 不 可 重 整 的 ,因此 就 需要 合适 的 正规 化 


机 制 消除 无 穷 大 。 这 里 使 用 三 动量 截断 ,动量 积 4 
上 限 为 4。 
NJL 模 型 只 有 夸克 自由 度 , 没 有 胶 子 自由 度 , 因 


此 也 就 没有 禁闭 效应 。 为 了 弥补 这 一 问题 ,此 处 考 
虑 了 Polyakov-loop, 即 BB 和 5B, 来 模拟 胶 子 的 部 分 效 
应 ,第 二 行 与 第 三 行 中 展示 了 其 与 夸克 的 相互 作用 。 
Tj o fU o Hp E TB EEA We HU HR, KG 
WA: 


73 Polyakov-loop 的 在 有 限 温 度 下 的 有 效 势 ， 
以 及 : 
J =( 


)(1 - 669 + 4(@ + d?) - 3(59)) 


(18) 
是 Polyakov-loop 到 它 的 迹 的 变换 的 雅克 比 行 
列 式 。 除 此 之 外 ,x 是 一 个 无 量 纲 的 参数 ,用 来 限制 
$$ 和 BB 的 值 在 0~1 之 间 。b,(7T) 是 温度 依赖 的 系数 ， 
其 形式 为 : 


T T 
b.(T) = a, + arp exp iy) (19) 


根据 文献 L50] , 模型 中 NJL 部 分 的 参数 可 以 通 
过 真空 的 性 质 来 确定 ,Polyakov-loop 部 分 的 参数 则 
可 以 通过 零 化 学 势 时 的 压强 密度 来 拟 合 ,这些 参数 
的 具体 值 展现 在 表 1 与 表 2 中 。 

在 平衡 态 系统 的 假设 下 S USER o, DL o IU 
B 是 自由 变量 ,其 值 可 以 通过 能 隙 方程 组 获得 : 


E 
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表 1 rPNJL 模 型 中 NJL 部 分 的 参数 
Table 1 Parameters for the NJL part in the rPNJL model 


ms/ MeV m,/ MeV A / MeV g4 


SpA g,/MeV* g,/MeV* 


5.5 183.468 637.720 2.914 


表 2 rPNJL 中 Polyakov-loop 部 分 参数 
Table 2 Parameters for the Polyakov-loop part in the 
rPNJL model 
T,/ MeV & a, a, b, b, K 
175 6.75 -98 0.26 0.805 7.555 0.1 


aQ aQ aQ 

do, ” ap 0G 

需要 注意 的 是 ,在 零 重 子 数 化 学 势 时 ,$B 和 FB 的 

值 相 等 ,但 是 在 有 限 化 学 势 时 则 不 相等 。 得 到 夸克 

BEER o, WA @ 和 5 在 给 定 温度 与 化 学 势 下 的 值 ,就 

可 以 通过 式 (5) 与 式 (6) 计 算 重 子 数 扰动 C,、C, 以 及 

kz。 在 零 化 学 势 时 ,利用 上 述 参数 计算 得 到 的 重子 

BPD xo? 随 温 度 的 变化 展示 在 图 3 中 ,其 中 77 = 

166 MeV ,为 IPNJL 给 出 的 零 化 学 势 下 的 相 变 温度 ， 

在 强 子 共 振 气 极限 下 xo? = 1, 在 自由 夸克 气 极限 时 
Ka” 约 0.068。 


0 (20) 


五 


m/m, 


TIT/ 


图 3 在 零 重 子 数 密度 下 利用 rPNJL 模型 计算 得 到 的 重子 数 

扰动 的 ka? 与 温度 的 关系 外 

Fig.3 Variation of the kurtosis of baryon number fluctuation 

kc? in the rPNJL model as a function of the temperature at zero 
baryon number density" 


图 3 中 虚线 为 IPNJL 模型 的 结果 ,不 同 的 点 取 
自 格 点 QCD Ait EO. BAB AA E 
其 中 ,7x = 166 MeV ArPNIL 模型 得 出 的 零 化 学 势 
下 的 相 变 温度 。 可 以 看 到 ,xo? 在 低温 下 趋向 于 1, 这 
正 是 HRG 模 型 给 出 的 低温 下 的 结果 ;而 在 高 温 下 则 
趋向 于 自由 夸克 气 的 结果 , 即 xa? 约 0.068。 整 体 上 
ITPNJL 模 型 计算 的 zc? 与 格 点 的 结果 比较 吻合 。 


4 ”结果 和 讨论 


为 了 确定 重子 数 扰动 xo* 随 对 撞 能 量 的 变化 关 

系 ,需要 先 确定 化 学 冻结 线 的 形式 。 此 处 选取 了 两 
条 化 学 冻结 线 的 形式 介 ED. mx: 
T (u) = 0.158 - 0.144? - 0.0444 


75.968 7 2.193 x 107 - 5.890 x 10°” 


这 个 形式 取 自 文献 [38] ,是 基于 了 RHIC 实验” 
数据 中 拟 合 的 ,fl 则 是 在 人 包 的 基础 上 做 了 微调 ,其 表 
达 式 为 : 

T (u) = 0.158 - 0.144? - 0.044 — 0.0lexp(— (u - 

0.067 )/0.05) 

式 中 :温度 7 和 重子 数 化 学 势 心 的 单位 都 是 GeV， 
下 同 。 

图 4 同时 展示 了 1PNJL 模 型 计算 得 到 的 相 图 
其 中 :CEP 的 位 置 为 (ys = 720 MeV , T? = 93 MeV), 
即 图 中 的 黑色 三 角形 处 ; 黑 线 为 一 级 相 变 线 ; 蓝 色 虚 
线 为 crossover; 浅 蓝 色 点 为 RHIC 的 BES-I 对 撞 结 
果 ; 圆 点 与 黑色 方块 分 别 为 RHIC BES-I5 5 以 及 低 
能 重 离子 对 撞 中 化 学 冻 出 的 温度 与 化 学 势 呈 ;同时 ， 
红色 与 蓝 色 虚 线 分 别 是 两 条 化 学 冻结 线 fl 和 亿 。 可 
以 看 出 ,fl 与 包 的 不 同 主要 体现 在 小 化 学 势 附近 ,也 
就 是 高 对 撞 能 量 下 :化 学 冻结 线 fl 从 相 边 界 的 “后 
方 ” 开 始 ,而 包 则 从 相 边 界 的 前方? 开始。 另外 ,为 
了 与 实验 结果 , 即 xo? 与 对 撞 能 量 的 关系 , 相 比较 ,还 
需要 重子 数 化 学 势 与 对 撞 能 量 的 关系 ,具体 形式 
如 下 : 


(/s)= 1.477 ji 
HEIN 1 + 0.343 /s Pn 


0 200 400 600 800 1 000 
Hz / MeV 


4 rPNJL EC TESERGAH EU" 
Fig.4 Chiral phase-transition line in the rPNJL model and 
freeze-out lines extracted from experiments"' 


ig 2X ,TPNJL VE 5 AY co? Ye a 45 7: UR GZ £1 FA £2. 
的 行为 展示 在 图 5 中, 分别 为 红色 虚线 与 蓝 色 虚线 ; 
绿色 虚线 为 NIL 计算 的 结果 5 ; 黑 点 为 STAR 的 净 
质子 数 测量 结果 ™”。 

首先 选用 化 学 冻结 线 f, 它 是 从 相 边界 的 “后 
方 开 始 , 可 以 看 到 , 沿 着 fl,rPNJL 模型 计算 的 重子 
数 扰动 xo? 随 着 对 撞 能 量 的 降低 ,首先 展现 出 了 一 个 
“下 四 ”结构 ,最 低 值 大 约 处 于 Vs = 20 GeV , 随 着 对 
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Vs /GeV 
图 5 ITPNJL 模型 计算 的 重子 数 xz? 随 对 撞 能 量 的 


变化 行为 加 
Fig.5 xo’ as a function of collision energy from the 
rPNJL model?" 


撞 能 量 继续 降低 ,在 WS = 6 GeV 处 出 现 一 个 峰 结 


化 学 冻结 
量 的 变化 关系 , 令 人 惊讶 的 是 ,计算 结果 基本 符合 实 


吉 线 方程 ,计算 了 重子 数 扰动 的 xo? 随 对 撞 能 


验 数 据 的 趋势 ,在 Vs = 7.7 GeV 处 产生 的 峰 结构 ， 


是 CEP 存 在 的 标志 。 考 虑 到 实验 观测 到 的 物理 量 


为 系统 化 学 冻 出 时 的 状态 ,而 此 时 系统 已 经 演化 到 


现 


平衡 态 ,因此 ,本 文中 基于 平衡 态 假 设 的 计算 结果 展 
出 与 实验 数据 较 好 的 符号 ,这 或 许 意味 着 在 理论 


icq ee a et 
合理 的 。 男 外 需要 强调 一 点 的 是 ,得 到 跟 实验 结 


果 相 符合 的 结果 ,最 主要 的 原因 是 以 下 用 点 :1 模型 
在 零 密 给 出 了 合理 的 相 变 温度 ,并 预言 了 比较 合理 
的 临界 终点 位 置 ,由 此 给 出 了 比较 合理 的 相 变 线 ;2) 


相 变 线 跟 化 学 冻结 线 的 位 置 关 系 对 最 终 xo? 随 对 撞 


能 量 的 变化 关系 非常 敏感 。 想 强调 的 一 点 是 


用 里 


ze Ko? [i 


Y , BE Jes o? 降低, 最 终 降 为 1 左右。 这 个 整体 的 行 
为 与 STAR 实验 测量 的 Vs > 7.7 GeV 的 结果 符合 得 
很 好 。 沿 着 第 二 条 化 学 冻结 线 亿 ,也 就 是 从 相 边 界 
的 前 方 出 发 ,kxo 只 有 在 Vs = 6Gev 产 生 一 个 峰 结 
构 , 而 没有 “下 四 ”结构 出 现 。 作 为 比较 ,文献 [10j 使 
用 的 NL 模型 计算 的 结果 与 实验 结果 差距 较 大 。 

此 可 以 看 出 ,重子 数 的 xo? 随 着 对 接 能 量 先 下 
降 后 上 升 的 行为 ,与 QCD 相 变 有 着 密切 的 关系 :在 
小 化 学 势 时 ( 即 高 能 对 撞 ) ,化 学 冻结 线 要 穿 过 


crossover 的 区 域 , 从 QGP 相 到 强 子 气相 ; 随 着 对 撞 
能 量 降低 ,化 学 冻结 线 逐 渐 靠 近 CEP JF PARC ER XC HE 


过 ,因此 kao’? 逐 渐 增 大 后 又 降低 ,形成 峰 结构 ,此 处 的 
峰 结 构 就 是 CEP 存 在 的 标志 。 

在 3 GeV 以 及 更 低 的 能 量 下 ,实验 结果 显示 了 
负 值 ,目前 通常 认为 这 是 由 于 低能 下 的 重子 数 守恒 
效应 导致 的 外。 这 个 结果 无 法 从 本 文中 使 用 的 
rPNJL 模型 得 出 ,从 图 $ 中 可 以 看 到 ,在 低 对 撞 能 量 
下 ,xc 趋向 于 1 而 不 是 负 值 , 这 说 明 TPNJL 模型 无 法 
描述 高 化 学 势 的 相 变 行为 和 相 结构 ,这 或 许 是 由 于 
TPNJL 模 型 中 只 有 夸克 自由 度 ,而 缺少 强 子 ,后 者 是 
强 子 气相 的 主要 自由 度 。 一 种 可 能 的 原因 是 , 强 子 
相 在 低温 有 限 化 学 势 处 存在 强 子 的 气 液 相 变 ,进而 
也 会 有 相应 临界 终点 ,在 低能 打靶 实验 中 ,产生 的 系 
统 的 化 学 势 较 大 而 温度 很 低 ,接近 强 子 气 液 相 变 的 
临界 终点 区 域 ,正如 QCD 手 征 相 变 CEP 附近 的 负 值 
区 域 类 似 ,此 时 也 有 可 能 出 现 xc? 为 负 的 结果 。 


5 结语 


由 于 守恒 荷 在 相 变 附近 会 出 现 剧烈 的 扰动 , 基 
此 ,在 重 离子 对 撞 中 可 以 通过 重子 数 扰动 来 寻找 
CEP。 使 用 3 味 的 PNJL 模型 ,其 在 零 化 学 势 时 重子 
数 扰动 与 格 点 结果 一 致 。 利 用 通过 实验 数据 拟 合 的 


[2:3 


E BE tec DICA E TE RE T8 B6 Bt 19.6 GeV 处 的 最 低 点 


JE H 


日 于 在 低 化 学 势 情况 下 化 学 冻结 线 高 于 相 变 线 决 


定 的 。 


作者 贡献 声明 


许 坤 完成 论文 初稿 的 写作 与 稿件 的 


修改 ;黄梅 提供 论文 撰写 思路 及 修改 稿件 。 
参考 文献 


1 


040005-6 


Fodor Z, Katz S D. A new method to study lattice QCD at 
finite temperature and chemical potential[J]. Physics 
Letters B, 2002, 534(1 - 4): 87 - 92. DOI: 10.1016/ 
80370-2693(02)01583-6. 

Ding H T, Karsch F, Mukherjee S. Thermodynamics of 
strong-interaction matter from lattice QCD[M]. Quark- 
Gluon Plasma 5. World Scientific, 2016: 1 - 65. DOI: 10. 
1142/9789814663717 0001. 

Schmidt C, Sharma S. The phase structure of QCD[J]. 
Journal of Physics G, 2017, 44(10): 104002. DOI: 10. 
1088/1361-6471/aa824a. 

Pisarski R D, Wilczek F. Remarks on the chiral phase 
transition in chromodynamics[J]. Physical Review D, 
1984, 29(2): 338 - 341. DOI: 10.1103/physrevd.29.338. 
Hatta Y, Ikeda T. Universality, the QCD critical and 
tricritical point, and the quark number susceptibility[J]. 
Physical Review D, 2003, 67: 014028. DOI: 10.1103/ 
physrevd.67.014028. 

Schwarz T M, Klevansky S P, Papp G. Phase diagram and 
bulk thermodynamical quantities in the Nambu-Jona- 
Lasinio model at finite temperature and density[J]. 
Physical Review C, 1999, 60(5): 055205. DOI: 10.1103/ 
physrevc.60.055205. 

Bowman E S, Kapusta J I. Critical points in the linear o 
model with quarks[J/OL]. Physical Review C, 2009, 79: 


202306.00141v1 


chinaXiv 


Ed 


F hij. QCD 临界 终点 与 


ChinaXiv 合 作 期 刊 


重子 数 扰动 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


015202. DOI: 10.1103/PhysRevC.79.015202. 

Mao H, Jin J S, Huang M. Phase diagram and 
thermodynamics of the Polyakov linear sigma model with 
three quark flavors[J]. Journal of Physics G: Nuclear and 
Particle Physics, 2010, 37(3): 035001. DOI: 10.1088/ 
0954-3899/37/3/035001. 

Schaefer B J, Wagner M. QCD critical region and higher 
moments for three-flavor models[J]. Physical Review D, 
2012, 85(3): 034027. DOI: 10.1103/physrevd.85.034027. 
Schaefer B J, Wagner M. Higher-order ratios of baryon 
number cumulants[J]. Central European Journal of 
Physics, 2012, 10(6): 1326 - 1329. DOI: 10.2478/s11534- 
012-0115-y. 

Qin S X, Chang L, Chen H, et al. Phase diagram and 
critical endpoint for  strongly-interacting quarks[J]. 
Physical Review Letters, 2011, 106(17): 172301. DOI: 
10.1103/PhysRevLett.106.172301. 

Luecker J, Fischer C S, Fister L, et al. Critical point and 
deconfinement from Dyson - Schwinger equations[C]. 
Proceedings of 8th International Workshop on Critical 
Point and Onset of Deconfinement (CPOD). 2013: 057. 
DOI: 10.22323/1.185.0057. 

Fu W J, Pawlowski J M, Rennecke F, et al. Baryon 
number fluctuations at finite temperature and density[J]. 
Physical Review D, 2016, 94(11): 116020. DOI: 10.1103/ 
PhysRevD.94.116020. 

Stephanov M A. Non-Gaussian fluctuations near the QCD 
critical point[J]. Physical Review Letters, 2009, 102(3): 
032301. DOI: 10.1103/PhysRevLett.102.032301. 
Stephanov M A. On the sign of kurtosis near the QCD 
critical point[J]. Physical Review Letters, 2011, 107(5): 
052301. DOI: 10.1103/PhysRevLett.107.052301. 
Asakawa M, Ejiri S, Kitazawa M. Third moments of 
conserved charges as probes of QCD phase structure[J]. 
Physical Review Letters, 2009, 103(26): 262301. DOI: 
10.1103/PhysRevLett.103.262301. 

Athanasiou C, Rajagopal K, Stephanov M. Using higher 
moments of fluctuations and their ratios in the search for 
the QCD critical point[J]. Physical Review D, 2010, 82 
(7): 074008. DOI: 10.1103/physrevd.82.074008. 

Zhu L L, Wang B, Wang M, et al. Energy and centrality 
dependence of light nuclei production in relativistic 
heavy-ion collisions[J]. Nuclear Science and Techniques, 
2022, 33(4): 45. DOI: 10.1007/s41365-022-01028-8. 
Wang M, Tao J Q, Zheng H, et al. Number-of-constituent- 
quark scaling of elliptic flow: a quantitative study[J]. 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


040005-7 


Nuclear Science and Techniques, 2022, 33(3): 37. DOL: 
10.1007/s41365-022-01019-9. 

Lan S W, Shi S S. Anisotropic flow in high baryon 
density region[J]. Nuclear Science and Techniques, 2022, 
33(2): 21. DOI: 10.1007/s41365-022-01006-0. 

Wang H, Chen J H. Anisotropy flows in Pb-Pb collisions 
at LHC energies from parton scatterings with heavy quark 
trigger[J]. Nuclear Science and Techniques, 2022, 33(2): 
15. DOI: 10.1007/s41365-022-00999-y. 

HSE, 黄 旭 光 . 相对 论 自 旋 流 体力 学 [J/OL]. 物理 学 报 ， 
2023. DOI: 10.7498/aps.72. 20230036. 

PU Shi, HUANG Xuguang. Relativistic spin 
hydrodynamics[J/OL]. Acta Physica Sinica, 2023. DOI: 
10.7498/aps.72.20230036. 


E 


Lin Z W, Zheng L. Further developments of a multi-phase 
transport model for relativistic nuclear collisions[J]. 
Nuclear Science and Techniques, 2021, 32(10): 113. DOI: 
10.1007/s41365-021-00944-5. 

Xu J F. Bulk viscosity of interacting magnetized strange 
quark matter[J]. Nuclear Science and Techniques, 2021, 
32(10): 111. DOI: 10.1007/s41365-021-00954-3. 

Su Y, Sun Y J, Zhang Y F, et al. Non-extensive statistical 
distributions of charmed meson production in Pb-Pb and 
pp(p) collisions[J]. Nuclear Science and Techniques, 
2021, 32(10): 108. DOI: 10.1007/s41365-021-00945-4. 
Wang H, Chen J H. Study on open charm hadron 
production and angular correlation in high-energy nuclear 
collisions[J]. Nuclear Science and Techniques, 2021, 32 
(1): 2. DOI: 10.1007/s41365-020-00839-x. 

Tang Z B, Zha W M, Zhang Y F. An experimental review 
of open heavy flavor and quarkonium production at RHIC 
[J]. Nuclear Science and Techniques, 2020, 31(8): 81. 
DOI: 10.1007/s41365-020-00785-8. 

Gao J H, Ma G L, Pu S, et al. Recent developments in 
chiral and spin polarization effects in heavy-ion collisions 
[J]. Nuclear Science and Techniques, 2020, 31(9): 90. 
DOI: 10.1007/s41365-020-00801-x. 

Liu Y C, Huang X G. Anomalous chiral transports and 
spin polarization in heavy-ion collisions[J]. Nuclear 
Science and Techniques, 2020, 31(6): 56. DOI: 10.1007/ 
s41365-020-00764-z. 

Wang F Q, Zhao J. Search for the chiral magnetic effect 
in heavy ion  collisions[J]. Nuclear Science and 
Techniques, 2018, 29(12): 179. DOI: 10.1007/s41365- 
018-0520-z. 

Hattori K, Huang X G. Novel quantum phenomena 


202306.00141v1 


chinaXiv 


E od 术 


ChinaXiv 合 作 期 刊 


2023, 46: 040005 


32 


33 


34 


35 


36 


37 


38 


39 


induced by strong magnetic fields in heavy-ion collisions 
[J]. Nuclear Science and Techniques, 2017, 28(2): 26. 
DOT: 10.1007/s41365-016-0178-3. 

孙 旭 , 周 晨 升 , 陈 金 辉 , 等 . 重 离子 磁 撞 中 QCD 物质 整体 
极 化 的 实验 测量 [J/OL]. 物理 学 报 , 2023. DOL: 10.7498/ 
aps.72.20222452. 

SUN Xu, ZHOU Chensheng, CHEN Jinhui, et al. 


Measurements of global polarization of QCD matter in 
heavy-ion collisions[J/OL]. Acta Physica Sinica, 2023. 
DOI: 10.7498/aps.72.20222452. 

高 建华 , 黄 旭 光 , 梁 作 堂 , 等 . 强 相 互 作 用 自 旋 轨道 耦合 
与 QGP 整 体 极 化 [MOL]. 物理 学 报 , 2023. DOI: 10.7498/ 
aps.72.20230102. 

GAO Jianhua, HUANG Xuguang, LIANG Zuotang, et al. 


Spin-orbital coupling in strong interaction and global spin 
polarization[J/OL]. Acta Physica Sinica, 2023. DOI: 
10.7498/aps.72.20230102. 

江 泽 方 , RRS, 余华 清 , 等 . RHIC HEX AutAu fled rp 
带电 粒子 直接 流 与 超 子 整体 极 化 的 计算 与 分 析 [JMOL]. 
物理 学 报 , 2023. DOI: 10.7498/aps.72.20222391. 

JIANG Zefang, WU Xiangyu YU Huaqing, et al. The 
the global 


direct flow of charged particles and 
polarization of hyperons in 200 AGeV Au-Au collisions 
at RHIC[J/OL]. Acta Physica Sinica, 2023. DOI: 10.7498/ 
aps.72.20222301. 

盛 欣 力 , 梁 作 堂 , 王 群 . 重 离子 碰撞 中 的 矢量 介子 自 旋 
HE 列 [WOL]. 物理 学报 ，2023. DOI: 10.7498/ 
aps.72.20230071. 


SHENG Xinli, LIANG Zuotang, WANG Qun. Global 


spin alignment of vector mesons in heavy ion collisions[J/ 
OL]. Acta Physica 2023. DOI: 10.7498/ 
aps.72.20230071. 

Adamczyk L, Adkins J K, Agakishiev G, et al. Energy 


Sinica, 


dependence of moments 
distributions at RHIC[J]. Physical Review Letters, 2014, 
112: 032302. DOI: 10.1103/PhysRevLett.112.032302. 

Aggarwal M M, Ahammed Z, Alakhverdyants A V, et al. 


of net-proton multiplicity 


Higher moments of net-proton multiplicity distributions at 
RHIC[J]. Physical Review Letters, 2010, 105: 022302. 
DOI: 10.1103/PhysRevLett.105.022302. 

Luo X F, Xu N. Search for the QCD critical point with 
fluctuations of conserved quantities in relativistic heavy- 
ion collisions at RHIC: an overview[J]. Nuclear Science 
and Techniques, 2017, 28(8): 112. DOI: 10.1007/s41365- 
017-0257-0. 

Abdallah M S, Aboona B E, Adam J, et al. Measurements 


40 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


040005-8 


of proton high order cumulants in "E 3 GeV AutAu 
collisions and implications for the QCD critical point[J]. 
Physical Review Letters, 2022, 128(20): 202303. DOI: 
10.1103/PhysRevLett.128.202303. 

Critelli R, Noronha J, Noronha-Hostler J, et al. Critical 
point in the phase diagram of primordial quark-gluon 
matter from black hole physics[J]. Physical Review D, 
2017, 96(9): 096026. DOI: 10.1103/physrevd.96.096026. 
Wu S J, Shen C, Song H C. Dynamically exploring the 
QCD matter at finite temperatures and densities: a short 
review[J]. Chinese Physics Letters, 2021, 38(8): 081201. 
DOI: 10.1088/0256-307x/38/8/081201. 

Shen C, Yan L. Recent development of hydrodynamic 
modeling in heavy-ion collisions[J]. Nuclear Science and 
Techniques, 2020, 31(12): 122. DOI: 10.1007/s41365- 
020-00829-z. 

Fang R H, Dong R D, Hou D F, et al. Thermodynamics of 
the system of massive Dirac fermions in a uniform 
magnetic field[J]. Chinese Physics Letters, 2021, 38(9): 
091201. DOI: 10.1088/0256-307x/38/9/091201. 

Peng H H, Zhang J J, Sheng X L, et al. Ideal spin 
hydrodynamics from the Wigner function approach[J]. 
Chinese Physics Letters, 2021, 38(11): 116701. DOI: 
10.1088/0256-307x/38/11/116701. 

Vovchenko V, Steinheimer J, Philipsen O, er al. Cluster 
expansion model for QCD baryon number fluctuations: 
no phase transition at u,/7 [J]. Physical Review D, 2018, 
97(11): 114030. DOI: 10.1103/physrevd.97.114030. 

Fu J H. Higher moments of net-proton multiplicity 
distributions in heavy ion collisions at chemical freeze-out 
[J]. Physics Letters B, 2013, 722(1 - 3): 144 - 150. DOI: 
10.1016/j.physletb.2013.04.018. 

Karsch F, Redlich K. Probing freeze-out conditions in 
heavy ion collisions with moments of charge fluctuations 
[J]. Physics Letters B, 2011, 695(1 - 4): 136 - 142. DOI: 
10.1016/j.physletb.2010.10.046. 

Fu J H. Higher moments of multiplicity fluctuations in a 
hadron-resonance gas with exact conservation laws[J]. 
Physical Review C, 2017, 96(3): 034905. DOI: 10.1103/ 
PhysRevC.96.034905. 

Bazavov A, Ding H T, Hegde P, et al. The QCD equation 
of state to O(u,°) from lattice QCD[J]. Physical Review 
D, 2017, 95(5): 054504. DOI: 10.1103/PhysRevD. 95. 
054504. 

Bhattacharyya A, Ghosh S K, Maity S, et al. 


Reparametrizing the Polyakov-Nambu-jona-lasinio model 


202306.00141v1 


chinaXiv 


许 坤 等 : 


QCD 临界 终点 与 


ChinaXiv& f'ERBTII 


liri 


ETA 


51 


52 


53 


54 


[J]. Physical Review D, 2017, 95(5): 054005. DOI: 10. 
1103/physrevd.95.054005. 

Li Z B, Xu K, Wang X Y, et al. The kurtosis of net baryon 
number fluctuations from a realistic Polyakov-Nambu- 
Jona-Lasinio model along the experimental freeze-out line 
[J]. The European Physical Journal C, 2019, 79(3): 245. 
DOI: 10.1140/epjc/s10052-019-6703-x. 

Ghosh S K, Mukherjee T K, Mustafa M G, et al. PNJL 
model with a Van der Monde term[J]. Physical Review D, 
2008, 77: 094024. DOI: 10.1103/PhysRevD.77.094024. 
Das S. Identified particle production and freeze-out 
properties in heavy-ion collisions at RHIC Beam Energy 
Scan program[J]. EPJ Web of Conferences, 2015, 90: 
08007. DOI: 10.1051/epjconf/20159008007. 

Adamczyk L, Adkins J K, Agakishiev G, et al. Bulk 
properties of the medium produced in relativistic heavy- 
ion collisions from the beam energy scan program[J]. 
Physical Review C, 2017, 96(4): 044904. DOI: 10.1103/ 
PhysRevC.96.044904. 


25 


56 


57 


58 


040005-9 


Kaczmarek O. Lattice QCD results on soft and hard 
probes of strongly interacting matter[J]. Nuclear Physics 
A, 2017, 967: 137 - 144. DOI: 10.1016/j. nuclphysa. 
2017.05.106. 

Begun V V, Vovchenko V, Gorenstein M I. Updates to the 
p*p and A+A chemical freeze-out lines from the new 
experimental data[J]. Journal of Physics: Conference 
Series, 2017, 779: 012080. DOI: 10.1088/1742-6596/779/ 
1/012080. 

Fan W K, Luo X F, Zong H S. Mapping the QCD phase 
diagram with susceptibilities of conserved charges within 
Nambu-Jona-Lasinio model[J]. International Journal of 
Modern Physics A, 2017, 32(11): 1750061. DOI: 10.1142/ 
s0217751x17500610. 

Braun-Munzinger P, Friman B, Redlich K, ef al. 
Relativistic nuclear collisions: establishing a non-critical 
baseline for fluctuation measurements[J]. Nuclear Physics 
A, 2021, 1008: 122141. DOI: 10.1016/j. nuclphysa. 2021. 
122141. 


